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Low-frequency  (<  10  Hz)  volcanic  earthquakes  originate  at  a  wide  range  of  depths  and  occur  before,  during,  and 
after  magmatic  eruptions.  The  characteristics  of  these  earthquakes  suggest  that  they  are  not  typical  tectonic  events. 
Physically  analogous  processes  occur  in  hydraulic  fracturing  of  rock  formations,  low-frequency  icequakes  in  temperate 
glaciers,  and  autoresonance  in  hydroelectric  power  stations.  We  propose  that  unsteady  fluid  flow  in  volcanic  conduits 
is  the  common  source  mechanism  of  low-frequency  volcanic  earthquakes  (tremor).  The  fluid  dynamic  source  mechan¬ 
ism  explains  low-frequency  earthquakes  of  arbitrary  duration,  magnitude,  and  depth  of  origin,  as  unsteady  flow  is  in¬ 
dependent  of  physical  properties  of  the  fluid  and  conduit.  Fluid  transients  occur  in  both  low-viscosity  gases  and 
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high-viscosity  liquids.  A  fluid  transient  analysis  can  be  formulated  as  generally  as  is  warranted  by  knowledge  of  the  com¬ 
position  and  physical  properties  of  the  fluid,  material  properties,  geometry  and  roughness  of  the  conduit,  and  boundary 
conditions.  To  demonstrate  the  analytical  potential  of  the  fluid  dynamic  theory,  we  consider  a  single-phase  fluid,  a  melt 
of  Mount  Hood  andesite  at  1250°C,  in  which  significant  pressure  and  velocity  variations  occur  only  in  the  longitudinal 
direction.  Further  simplification  of  the  conservation  of  mass  and  momentum  equations  presents  an  eigenvalue  problem 
that  is  solved  to  determine  the  natural  frequencies  and  associated  damping  of  flow  and  pressure  oscillations.  For  a  simple 
constant  pressure  reservoir-conduit-orifice  fluid  system,  a  change  in  orifice  size  can  cause  the  source  signal  duration  due 
to  a  single  disturbance  to  change  from  several  seconds  to  in  excess  of  an  hour.  Fluid  kinematic  viscosity  (0.412  m2/s)  is 
both  included  and  neglected  in  the  analysis  to  illustrate  its  effect  upon  system  damping.  Tremor  magnitude  is  related 
to  the  magnitude  of  the  conduit  wall  motion  occurring  in  response  to  oscillating  fluid  pressures.  The  entire  volcanic 
fluid  system  does  not  generally  experience  a  transient  flow  condition  at  a  given  time.  Closed  end  and  constant  pressure 
fluid  system  boundaries  reflect  pressure  waves  and  isolate  components  of  the  fluid  system.  Total  fluid  movement  may 
not  be  accurately  monitored  by  tremor  data  alone,  since  steady  state  and  slowly  changing  flows  are  aseismic. 
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FLUID  DYNAMIC  ANALYSIS  OF 
VOLCANIC  TREMOR 

M.C.  Ferrick,  A.  Qamar  and  W.F.  St.  Lawrence 


INTRODUCTION 

The  characteristics  of  volcanic  earthquakes  sug¬ 
gest  that  they  are  not  typical  tectonic  events.  Be¬ 
ginning  with  Omori  (1914),  many  authors  (e.g. 
Kubotera  1974,  Minakami  1974a,  Tanaka  1979)  have 
described  and  classified  volcanic  earthquakes  accord¬ 
ing  to  frequency  content  and  duration  of  the  seismic 
radiation.  In  general,  investigators  refer  to  three 
types  of  seismic  signals  produced  by  volcanic  earth¬ 
quakes  that  are  not  directly  associated  with  erup¬ 
tions:  type  A  and  B  events  and  harmonic  tremor. 
High-frequency  (>  10-Hz)  type  A  events  have  seismic 
signatures  like  those  of  tectonic  earthquakes.  Type 


B  events  yield  signals  having  a  low-frequency  con¬ 
tent  in  the  range  of  0  to  5  Hz.  Sustained  low-frequency 
(0  to  10  Hz)  ground  motions  lasting  a  few  minutes  to 
several  hours  are  termed  harmonic  tremor.  Figure  1 
includes  sample  signatures  of  a  type  B  event  and  a 
portion  of  a  harmonic  tremor. 

Aki  et  al.  (1 977)  have  noted  that  low-frequency 
volcanic  earthquakes  have  many  characteristics  that 
are  similar  to  those  observed  in  seismic  signals  emitted 
from  the  hydraulic  fracturing  experiments  of  Potter 
and  Dennis  (1974).  Characteristics  of  volcanic  tremor 
(Aki  et  al.  1977,  Aki  and  Koyanagi  1981,  Chouet 
1981,  Seidl  et  al.  1981)  include  a  stationary  seismic 
source,  a  peaked  spectrum  and  a  dominant  frequency 
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a.  Type  B  volcanic  earthquake  (the  large  signal). 


lO^im 


ISt 


b.  Two  one-minute  segments  of  a  harmonic  tremor  recorded  five  minutes  apart.  This  tremor 
lasted  for  about  20  minutes.  The  ground  motion  scale  was  calculated  at  a  frequency  of  1.5 
Hz.  In  this  figure  a  beating  pattern  can  be  noted. 


Figure  1,  Seismograms  recorded  at  the  Tlmberllne  Parking  Lot  on  Mount  St.  Helens  on 
4  April  1980. 


that  can  change  during  the  course  of  eruptive  activ¬ 
ity. 

The  origin  of  low-frequency  volcanic  earthquakes 
is  controversial.  Some  theories  of  the  origin  of  the 
tremors  are  free  oscillations  of  a  magma  chamber 
(Kubotera  1974,  Sassa  1935,  Shima  1958),  oscilla¬ 
tions  of  gas  in  a  volcanic  vent  (Steinberg  and  Stein¬ 
berg  1975),  oscillation  of  volcanic  layers  from  mov¬ 
ing  magma  (Omer  1950)  and  random  opening  of 
tensile  cracks  from  excess  fluid  pressure  (Aki  et  al. 
1977,  Aki  and  Koyanagi  1981,  Chouet  1981).  All 
of  these  theories  propose  source  effects  as  the  cause 
of  tremor.  The  observation  of  predominant  low- 
frequency  ground  motion,  even  at  epicentral  dis¬ 
tances  of  only  a  few  kilometers,  is  evidence  that 
the  motion  is  produced  at  the  source  rather  than 
being  an  effect  of  wave  propagation  through  a  heter¬ 
ogeneous  medium.  Seidl  etal.  (1981)  noted  the 
stability  of  low-frequency  peaks  obtained  at  the 
recording  stations  near  Mount  Etna  and  concluded 
that  the  origin  was  source  related. 

Volcanic  earthquakes  are  known  to  occur  at  ba¬ 
saltic,  andesitic  and  dacitic  volcanoes  as  forerunners 
of  eruptions  (Minakami  1974b).  Changes  in  seismic 
activity  have,  for  example,  been  one  of  the  consis¬ 
tent  precursors  of  volcanic  eruptions  at  Mount  St 
Helens,  Washington.  Over  a  thousand  volcanic  earth¬ 
quakes  exceeding  magnitude  3,  and  at  least  18  epi¬ 
sodes  of  harmonic  tremor,  occurred  before  the  major 
eruption  of  18  May  1980.  The  most  violent  tremors 
occurred  several  hours  after  the  initial  explosion,  and 
since  that  time  the  earthquakes  have  been  much 
smaller.  A  marked  increase  in  seismicity  and  one  or 
more  episodes  of  harmonic  tremor  have  preceded 
and  accompanied  most  subsequent  eruptions  (Decker 
and  Decker  1981,  Lipman  and  Mullineaux  1981). 

In  this  analysis  we  seek  an  understanding  of  the 
source  mechanism  of  both  low-frequency  harmonic 
tremor  and  volcanic  (type  B)  earthquakes.  Such  an 
understanding  allows  a  framework  for  interpretation 
of  many  apparently  diverse  observations.  Physical 
analogies  identify  unsteady  fluid  dynamics  as  a  likely 
source  mechanism  of  tremor.  Analyses  of  tremor  to 
date  have  typically  relied  upon  hypothesized  fluid 
motions.  Unlike  previous  investigation,  we  consider 
directly  the  dynamics  of  the  fluid  system.  Equations 
describing  the  fluid  vibration  are  developed  from  the 
physical  laws  of  conservation  of  mass  and  momentum. 
Fluid  properties  characteristic  of  Mount  St  Helens 
magmas  are  used  to  demonstrate  die  plausibility 
and  analytical  potential  of  the  fluid  dynamic  source 
mechanism. 


CONCEPT 

Steady  fluid  flow,  by  definition,  occurs  with  time- 
invariant  fluid  pressure  at  all  points  in  the  system.  A 
fluid  system  at  rest,  for  example,  is  at  a  steady  (zero) 
flow.  Excitations  that  occur  slowly  relative  to  the 
fundamental  period  of  the  fluid  system  induce  gradual 
(aseismic)  pressure  variations.  When  a  fluid  system 
initially  at  steady  state  is  abruptly  disturbed,  a  per¬ 
iod  of  transient  flow  results.  All  fluid  systems,  inde¬ 
pendent  of  conduit  geometry  and  fluid  properties, 
experience  fluid  transients  when  perturbed.  Appli¬ 
cation  of  these  concepts  to  volcanoes  requires  only 
that  we  presume  the  existence  of  a  fluid  system. 

We  propose  that  low-frequency  volcanic  events 
of  various  durations  (volcanic  tremor)  have  unsteady 
fluid  flow  as  a  common  origin.  The  initial  excitation 
to  which  the  fluid  system  responds  could  be  provided 
by  a  number  of  mechanisms.  What  is  required  is  an 
abrupt  flow  input,  an  abrupt  outflow,  or  some  other 
perturbation  of  an  existing  steady  condition.  A  likely 
mechanism,  important  during  pre-eruptive  periods, 
is  the  fracture  of  a  conduit  and  the  subsequent  fluid 
response  of  a  system  initially  at  rest.  Once  initiated, 
flow  and  pressure  oscillations  continue  until  damped 
as  a  result  of  fluid  friction,  generation  of  seismic  ra¬ 
diation  and  outflow  of  oscillatory  energy  from  the 
system.  The  pressure  oscillations  in  the  fluid  drive  an 
oscillatory  displacement  of  the  conduit  wall  and 
generate  elastic  waves  in  the  surrounding  rock. 

This  source  mechanism  is  suggested  by  a  number 
of  physical  analogies.  On  17  June  1976,  the  Edward 
Hyatt  Power  Plant  at  Oroville  Dam,  Oroville,  Cali¬ 
fornia,  experienced  severe  vibrations  in  penstock  no. 

1  (Ehrhart  1979).  A  transient  flow  condition  in  the 
penstock  with  severe  pressure  oscillations  at  a  natural 
frequency  of  the  fluid  system  was  produced  and  main¬ 
tained  (forced)  by  a  malfunction  of  the  turbine  shut¬ 
off  valve.  The  low-frequency  seismic  signals  emanating 
from  the  penstock  were  recorded  by  the  network  of 
seismometers  in  the  area.  Seismic  activity  ceased 
when  the  vibration  was  arrested.  McNutt  (1982)  re¬ 
ported  similar  observations  of  low-frequency  seis¬ 
micity  within  50  km  of  Tarbela  Dam,  Pakistan.  He 
compared  the  seismic  spectra  near  Tarbela  Dam  with 
the  tremor  spectra  at  six  circum-Pacific  stratovolcanoes 
and  found  them  to  be  similar,  each  exhibiting  narrow, 
evenly  spaced  peaks.  An  important  observation 
addressed  analytically  in  this  paper  was  that  the  re¬ 
corded  spectra  near  Tarbela  Dam  changed  with  the 
size  of  the  outflow  opening  of  the  tunnel. 

St  Lawrence  and  Qamar  (1979a,  1979b)  observed 
that  low  frequency  glacial  icequakes  and  volcanic 


2 


b.  After  hydraulic  fracture  treatment. 

Figure  2.  Time  variation  of  fluid  pressure  near  the  surface  of  a 
well  of  1590-m  depth  and  0. 062-m  diameter.  The  fluid  transient 
was  induced  by  rapid  fluid  withdrawal  at  the  surface  in  each  in¬ 
stance.  The  frequency  of  the  oscillation  has  decreased  by  more 
than  a  factor  of  two  between  these  pulse  tests.  Static  pressure 
Is  greatly  exceeded  during  each  transient  flow  condition. 
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earthquakes  produce  seismic  signals  that  cannot  be 
differentiated  from  one  another.  A  common  link 
between  these  two  physical  systems  is  that  each  in¬ 
cludes  fluid-filled  conduits. 

Anderson  and  Stahl  (1967)  conducted  field  tests 
to  investigate  the  process  of  hydraulic  fracturing. 
Using  formation  crude  oil  as  the  fracture  fluid,  they 
observed  fluid  pressure  oscillations  in  the  wellbore. 
Subsurface  fluid  pressure  transients  generated  during 
hydraulic  fracturing  of  the  formation  exhibited  the 
same  characteristics  both  deep  in  the  well  and  near 
the  surface.  Fracturing  of  the  formation  caused  the 
period  of  the  oscillations  to  increase. 

Recently  we  conducted  experiments  in  conjunc¬ 
tion  with  a  hydraulic  fracture  treatment  of  a  rock 
formation  near  the  bottom  of  a  1 590-m  well  of 
0,062-m  diameter.  The  fluid  that  was  used  had  a 
relatively  high  viscosity  of  about  80  centipoise,  which 
together  with  the  small  diameter  conduit  gave  a 


highly  damped  response.  Testing  included  pre-fracture 
and  post-fracture  "pulse”  tests  in  which  a  small  volume 
of  fluid  was  abruptly  released,  by  valving  at  the  surface, 
from  the  pressurized  fluid  system  which  was  initially 
at  rest.  Fluid  pressure  variations  with  time  at  a  point 
near  the  surface,  as  a  result  of  the  fluid  release,  are 
presented  in  Figure  2  for  pre-  and  post-fracture  tests. 

The  pre-fracture  fluid  system  consisted  of  a  con¬ 
duit  which  was  closed  at  both  ends  so  that  zero  flow 
boundaries  were  present.  The  pressure-time  trace 
for  this  case  exhibits  a  dominant  frequency  of  0.44 
Hz.  Increasing  pressure  near  the  surface  is  accompan¬ 
ied  by  fluid  flow  toward  the  surface  and  conversely 
for  decreasing  pressure. 

After  the  fracture  treatment  of  the  formation,  a 
large  volume  of  fluid  remained  in  the  system  of  frac¬ 
tures  which  intersected  the  well  near  its  base.  The 
post-fracture  pulse  test  yielded  a  dominant  frequency 
of  0. 1 7  Hz.  The  downhole  boundary  had  changed 
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from  a  zero  flow  boundary  to  a  constant  pressure 
“reservoir”  boundary.  This  change  results  in  a 
doubling  of  the  fundamental  period  of  the  fluid  sys¬ 
tem.  The  post-fracture  period  was  greater  than  twice 
the  pre-fracture  period,  indicating  that  the  conduit 
length  has  increased.  The  post-fracture  trace  exhibts 
line  packing,  a  characteristic  of  high  friction  fluid 
systems.  Line  packing  occurs  when  flow  is  not 
totally  stopped  by  the  passage  of  the  compression 
wave;  the  pressure  near  the  closed  end  of  the  pipe 
continues  to  rise  for  a  period  of  time  following  the 
abrupt  pressure  increase.  In  each  pulse  test,  fluid 
pressures  recorded  near  the  surface  during  the  tran¬ 
sient  flow  condition  greatly  exceeded  the  existing 
static  pressure. 

Seidl  et  al.  (1981),  following  Sl  Lawrence  and 
Qamar  (1979a,  1979b),  have  applied  basic  concepts 
of  the  fluid  transient  theory  to  the  analysis  of  tremor 
at  Mount  Etna.  Analysis  of  seismometer  records 
near  the  volcano  revealed  that  eruptions  can  be  dis¬ 
tinguished  according  to  different  eigenfrequencies 
of  the  fissure  and  conduit  network.  Positions  of 
active  dykes  and  conduits  were  traced  following  the 
amplitude  distributions  associated  with  pronounced 
spectral  lines  in  the  tremor  spectra.  In  addition,  it 
was  found  that  harmonic  tremor  and  volcanic  (type 
B)  earthquakes  have  similar  spectra. 

LeGuern  et  al.  (1982)  observed  several  periodic 
swarms  of  between  10  and  14  emissions  of  gas  from 
a  vent  on  Mount  Etna.  Each  swarm  had  a  period  of 
1.8  s,  decreased  in  intensity  with  each  succeeding 
emission  in  the  swarm,  and  was  typically  followed 
by  a  pause  of  a  few  minutes.  Flow  velocity,  dynamic 
pressure  and  temperature  measured  at  the  vent  re¬ 
vealed  that  the  phenomenon  was  due  to  resonance 
of  the  volcanic  gas  in  the  conduit. 


BACKGROUND 

Pressure  and  flow  oscillations  occur  in  fluid  sys¬ 
tems  at  a  frequency  that  depends  upon  both  the 
geometry  of  the  container  and  the  location  of  the 
excitation  (Muto  and  Kanei  1980,  Zielke  and  Hack 
1972).  The  excitation  need  only  be  temporary  to 
induce  a  fluid  oscillation. 

Analysis  of  the  fluid  system  of  a  specific  volcano 
may  require  that  a  multiphase  fluid  (liquid,  gas,  solid 
in  suspension)  having  significant  temperature  varia¬ 
tions  be  considered.  The  general  theory  is  complex 
and  full  exposition  is  unwarranted  in  this  discussion. 
Therefore,  we  will  consider  a  single  phase  fluid 
(liquid)  in  which  temperature  variability  does  not 
significantly  affect  the  flow.  The  flow  of  a  gas  could 
also  be  considered  wilh  relatively  minor  modifications 
of  this  development  (Wylie  and  Streeter  1 978).  The 


nature  of  a  fluid  transient  in  a  closed  conduit  is  then 
dependent  upon  fluid  density,  viscosity  and  pressure 
wave  speed,  length  and  shape  of  the  conduit,  bound¬ 
ary  conditions  and  the  nature  of  the  disturbance. 

The  physical  laws  that  describe  the  flow  of  a 
liquid  in  a  conduit  are  (Wylie  and  Streeter  1978) 
the  conservation  of  momentum 
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where 


H  =  the  pressure  head  in  meters, 

<?  =  the  volumetric  flow  rate  (m3/s), 

A  =  the  cross-sectional  area  of  the  conduit 

(m2), 

F(Q)  =  a  frictional  loss  function  which  is  de¬ 
pendent  upon  fluid  viscosity,  conduit 
size,  shape,  and  roughness, 
a  =  the  angle  of  the  conduit  with  the  hori¬ 
zontal  , 

x  =  linear  distance  along  the  conduit, 
g  =  the  acceleration  due  to  gravity, 
a  =  the  pressure  wave  speed  in  the  field, 
t  =  time. 


These  quasi-linear  hyperbolic  equations  together  with 
appropriate  boundary  and  initial  conditions  can  be 
solved  numerically.  The  governing  equations  as  stated 
require  the  assumption  that  the  conduit  material  be¬ 
haves  elastically.  The  generation  of  seismic  radiation 
is  accomplished  without  energy  loss  by  the  fluid  sys¬ 
tem.  If  the  elastic  assumption  cannot  be  made,  a 
stress-strain  model  describing  conduit  material  be¬ 
havior  is  required  and  the  corresponding  fluid  dynamic 
model  would  exhibit  increased  damping. 

For  hyperbolic  equations,  H  and  Q  information 
propagates  in  the  x-t  plane  along  characteristics.  The 
slope  of  a  characteristic  in  closed  conduit  flow  appli¬ 
cations  is  dominated  by  die  pressure  wave  speed  in 
the  fluid: 


dx  Q  ^  . 

ST’a1"'1"- 


(3) 


Both  flow  and  pressure  at  all  points  in  the  fluid  re¬ 
spond  rapidly  to  a  disturbance  in  the  system. 

Equations  1  and  2  are  typically  the  basis  of  anal¬ 
yses  of  transients  in  fluid  systems  for  which  the  sys¬ 
tem  geometry,  boundary  conditions  and  fluid  proper¬ 
ties  are  known.  Given  sufficient  data  at  a  specific 
volcano,  the  behavior  of  the  fluid  transient  source 
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could  be  modeled  in  a  precise  way  by  using  these  or 
more  general  equations.  Typically,  however,  the 
configuration  of  the  fluid  system  is  poorly  known. 

We  are  then  faced  with  th$  inverse  problem  of  finding 
the  fluid  system  geometry  and  the  precise  nature  of 
the  frictional  loss  function.  The  solution  of  this 
problem  is  not  unique.  Volcanoes  with  different 
magmas  and  conduit  system  geometries  could  gen¬ 
erate  similar  low-frequency  seismic  signals. 

As  steady  fluid  flow  is  aseismic,  it  usually  is  not 
possible  to  quantify  the  total  volume  of  magma  flow 
using  seismic  data  generated  during  episodes  of  vol¬ 
canic  tremor.  If  low-frequency  pressure  oscillations 
at  the  source  can  be  determined  with  reasonable 
confidence,  however,  fluid  dynamic  modeling  tools 
can  be  applied  to  investigate  both  the  rate  of  oscil¬ 
latory  fluid  flow  and  the  physical  characteristics 
of  the  fluid  system. 


=  -L  +  S*tL  &L  ( 7 ) 

a*2  "  a2  bt2  a2  dt 

With  the  assumption  that/»'=  X(x)T(t),  the  method 
of  separation  of  variables  can  be  applied  to  eq  7 
yielding 


where  y2  is  the  separation  constant  The  particular 
solution  for  an  oscillatory  pressure  head  can  be  found 
from  eq  8  as 

h'  =  est  (c,*?*  +  C2ei'jrj 

y2  =  -Z-(s+gAR)  (9) 

a 1 


LINEAR  THEORY 

The  basic  features  of  a  fluid  system  response  to 
a  disturbance,  including  the  natural  frequencies  of 
the  vibration  and  associated  damping  of  the  oscilla¬ 
tions,  can  be  conveniently  determined  with  a  free 
vibration  analysis  (Wylie  and  Streeter  1978,  Zielke 
and  Hack  1972).  The  nonlinear  terms,  Q/A 2  bQ/bx 
and  Q/A  bH/bx  in  eq  1  and  2  respectively,  are  disre¬ 
garded  in  this  analysis  as  fluid  velocity  and  the  spa¬ 
tial  gradients  of  head  and  discharge  are  typically 
small.  The  frictional  loss  term,  F(Q),  which  is  non¬ 
linear  for  turbulent  flow,  must  be  linearized.  Also, 
for  cases  where  the  volumetric  flow  rate  is  not  ex¬ 
tremely  large,  the  sin  a  term  in  the  conservation  of 
mass  equat,v  is  small  and  can  be  neglected.  Finally, 
the  dependent  variables,  H  and  <?,  are  decomposed 
into  mean  (//,  Q)  and  fluctuating  (h’ ,q')  compon¬ 
ents. 

The  linearized  equations  for  oscillatory  flow  are 
then 


where  Cj  and  C2  are  constants. 

The  complex  frequency  s  is  a  constant  composed 
of  a  real  part  a  which  quantifies  system  damping,  and 
an  imaginary  part  w  the  circular  frequency  of  the 
oscillation.  Substituting  eq  9  into  eq  4  and  5  yields 
the  solution  for  q’ 

4 '  =  S~  *st  (c,  -  C2e~r*j  (1 0) 

If  the  upstream  end  of  the  conduit  is  identified  with 
subscript  U,  the  constants  C,  and  C2  are  determined 
as 

C,  =  i(//u-ze<?u) 

C2  =  \  (Hu  +ZC  Qu) 


7  -a—  1 

c~  gA  s 


(ID 


The  equations  for  h' and  q  can  then  be  written 


1 Rq,=0 
bx  gA  bt 

(4) 

h'  = 

§2.'+2d^'  =  0 

ax  a2  bt 

(5) 

<7- 

H, 


U  . 

-  sin 
c 


(12) 


where  R  is  the  linearized  resistance  per  unit  length. 
These  equations  can  be  cast  in  the  form  of  the  tele¬ 
graph  (wave)  equation  for  both  h‘  and  q‘: 


*hLsi-*?s:+sA**s: 

bx2  a2  bt2  a2  dt 


(6) 


The  expressions  for  H(x)  and  Q(x)  are  termed  trans¬ 
fer  equations  for  head  and  discharge  and  can  be  used 
to  relate  H  and  Q  at  each  end  of  a  conduit 

To  analyze  a  given  fluid  system,  transfer  equations 
for  all  lines  and  additional  equations  for  terminal 
conditions,  conduit  branch  points  and  other  compon¬ 
ents  are  assembled.  The  resulting  system  of  transfer 
equations  is  then  written  as 
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(13) 


Msx  =  0 

where  Ms  is  the  overall  system  matrix,  and  x  is  a 
vector  of  pressures  and  flows.  The  values  of  s  for 
which  the  matrix  equation  has  a  nontrivial  solution 
are  those  for  which  the  determinant  of  Mi  vanishes. 
The  residual  time  response  of  a  fluid  system  that 
has  undergone  an  excitation  is  known  when  the  com¬ 
plex  frequency  is  found  by  solution  of  this  eigen¬ 
value  problem  (Zielke  and  Hack  1972).  Because 
the  fluid  system  is  continuous,  an  infinite  number 
of  complex  frequencies  can  be  found.  Practically, 
however,  only  the  lowest  few  harmonics  are  ob¬ 
served. 

A  consequence  of  the  linear  theory  is  that  un¬ 
sustained  disturbances  in  fluid  systems  decay  ex¬ 
ponentially  with  time  ( eq  12).  Envelopes  of  dimen¬ 
sionless  amplitude  decay  are  plotted  in  Figure  3  as 
a  function  of  time.  Depending  upon  the  magnitude 
of  o,  oscillations  can  persist  for  a  few  seconds  or  in 
excess  of  an  hour  as  the  result  of  a  single  disturbance. 

Two  parameters  of  this  analysis  are  important  for 
characterizing  the  response  of  the  fluid  system  of 
the  volcano.  The  characteristic  impedance,  Zc,  is  a 
function  of  fluid  and  conduit  properties.  In  large 
diameter  pipes,  the  value  of  the  characteristic  impe¬ 
dance  approaches  that  for  frictionless  sytems  (eq  9 
and  11): 


Tim*  (*) 


0. 


The  fracture  of  the  pipe  wall  becomes  an  orifice 
through  which  fluid  flows  out  of  the  system.  The 
orifice  impedance  Zor  is  defined  as 

*or=^  (15) 

Q 

where  H  is  the  mean  pressure  head  and  Q  is  the 
mean  flow  just  upstream  of  the  orifice.  A  relatively 
large  opening  in  the  conduit  through  which  fluid 
passes  easily  with  only  a  small  loss  of  pressure  head 
has  a  small  orifice  impedance.  A  small  conduit  open¬ 
ing  which  does  not  readily  pass  fluid  has  a  large  ori¬ 
fice  impedance. 

The  frequency  and  damping  characteristics  of 
the  fluid  system  are  greatly  influenced  by  the  rela¬ 
tive  magnitudes  of  the  characteristic  and  orifice  im¬ 
pedances.  For  a  frictionless  reservoir-conduit-orifice 
fluid  system,  Hu  =0  and  R  =  0.  Identifying  the 
downstream  end  of  the  conduit,  near  the  orifice, 
with  subscript  D, 


Equations  12,  14-15  then  yield 

Zor  cosh  yf  +  Zc  sinh  yl  =  0  (1 7) 

where  /  is  the  conduit  length.  Solution  of  the  im¬ 
aginary  part  of  eq  17  gives  the  circular  frequencies 


Tim*  (min) 


b. 


Figure  3.  Envelopes  of  dimensionless  amplitude  decay  of  an  oscillation  for  a  range  of  system  damping  values  (a).  H0 
Is  the  original  or  maximum  amplitude  of  the  oscillation.  Note  the  change  of  scale  of  the  time  axis  between  Figures  3a 

and  3b. 
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Figure  4.  Single  conduit  components  of  a  fluid  system.  Intersections 
with  a  significantly  larger  conduit,  a  large  volume  fracture  zone  or  a 
dosed  end  are  boundary  conditions  that  Isolate  components  of  the 
fluid  system. 


of  the  system  and  the  real  part  of  eq  17  is  solved  to 
obtain  system  damping.  When  Zor  >  Zc,  the  fluid 
system  responds  as  if  the  conduit  were  a  closed  pipe 
and  the  frequencies  are  the  odd  harmonics  of  the 
conduit: 


Co  =  2irf  = 


nna  ,  _ 

2i_  i,.jfj,... 


o  = 


^Wor  ’ 


(18) 


If  Zc  >  Zor,  the  system  responds  as  if  the  conduit 
were  open-ended  and  the  frequencies  are  the  even 
harmonics  of  the  conduit: 


_  ,  nna  .  . , 

03  =  2nf=  ~2£,  n- 2,4,6... 


o  = 


Zc-Zq, 
Zc+Zor  ‘ 


(19) 


For  both  cases,  system  damping  increases  as  the  mag¬ 
nitudes  of  the  characteristic  and  orifice  impedance 
approach  each  other.  Less  pressure  wave  reflection 
occurs  at  the  orifice  and  oscillatory  energy  passing 
out  of  the  system  increases.  As  the  ratio  of  the  mag¬ 
nitudes  departs  from  one,  more  energy  is  reflected 
and  system  damping  decreases. 

The  dominant  range  of  tremor  frequencies  ob¬ 
served  at  a  number  of  volcanoes  is  1  to  2  Hz.  The 
natural  frequencies  of  a  fluid  transient  are  a  function 
of  the  pressure  wave  speed  in  the  fluid  and  of  conduit 
length  and  boundary  conditions  (eq  18).  Much  lower 
frequencies  would  be  expected  if  die  complete  fluid 
system  within  the  volcano  wen  experiencing  a  tran¬ 


sient  condition  at  a  given  time.  A  segment  of  the 
system  isfluid-dynamically  isolated  from  the  re¬ 
mainder  of  the  fluid  system  when  die  boundary  con¬ 
dition  provides  total  pressure  wave  reflection.  This 
can  occur  physically  by  conduit  pinchout  given  a 
closed  end  boundary  condition,  or  by  intersection 
with  a  significantly  larger  conduit  or  a  large  volume 
fracture  zone  creating  a  reservoir  or  constant  pressure 
boundary  condition  (Fig.  4).  The  complete  volcanic 
fluid  system  is  visualized  as  a  tree-like  system  of  these 
components.  Any  component  of  the  system  can  ex¬ 
perience  a  transient  flow-induced  vibration. 


APPLICATION  OF  LINEAR  THEORY: 

MOUNT  ST.  HELENS 

To  demonstrate  the  potential  insights  that  can  be 
obtained  from  the  fluid  dynamic  theory,  we  will  now 
make  an  assessment  of  some  features  of  volcanic 
tremor  at  Mount  St  Helens.  The  fluid  system  seg¬ 
ment  to  be  analyzed  is  assumed  to  be  composed  of 
a  relatively  large  body  of  magma  connected  to  a 
smaller  conduit  which  terminates  in  a  closed  end. 

The  larger  magma  body  could  be  a  chamber,  a  large 
conduit,  or  an  extensive  zone  of  fractures  containing 
a  volume  of  magma  that  is  large  relative  to  the  vol¬ 
ume  contained  in  the  conduit  itself.  The  smaller 
conduit  is  assumed  to  be  nearly  vertical,  with  an 
approximately  circular,  prismatic  cross  section  and  a 
diameter  of  100  m.  The  effect  of  choosing  a  relatively 
large  diameter  conduit  is  to  decrease  system  damping 
due  to  fluid  friction.  The  same  fluid  in  a  smaller 
conduit  would  exhibit  greater  viscous  damping.  The 
cylindrical  conduit  shape  is  chosen  to  take  full  ad¬ 
vantage  of  existing  fluid  transient  theory.  A  relatively 
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short  conduit  length  is  indicated  by  the  observed 
tremor  frequencies  (Fig.  1).  An  "order  of  magni¬ 
tude”  conduit  length  of  1000  m  is  chosen  for  this 
analysis. 

Murase  and  McBirney  (1973)  experimentally  de¬ 
termined  the  fluid  density  p,  kinematic  viscosity  u, 
bulk  modulus  K,  and  pressure  wave  speed  of  melts 
of  Mount  Hood  andesite  composed  of  6096  Si02. 

The  melt  was  totally  fluid  at  temperatures  greater 
than  1 200°C.  For  an  assumed  temperature  of  1 250°C, 
the  parameter  values  were  given  as  p  =  2430  kg/m3, 
v  =  0.41 2  m2/s,  K  =  1.78  x  1 0*°  Pa,  and  the  wave 
speed  in  an  unbounded  fluid,  a  -  2700  m/s.  it  is 
reasonable  to  expect  that  the  degassed  melt  had  a 
much  higher  viscosity  than  in-situ  magma  at  the  same 
temperature.  The  decreased  viscous  damping  due  to 
the  lower  in-situ  viscosity  is  neglected  in  this  analysis. 

In  closed  conduit  flow,  conduit  properties  affect 
the  pressure  wave  speed  of  the  fluid.  If  we  again 
assume  that  the  conduit  is  in  an  elastic  medium  and 
is  constrained  from  longitudinal  movement,  the  wave 
speed  can  be  expressed  as  (Wylie  and  Streeter  1978): 


rj2  = 


Kip 


(20) 


where  p  and  £R  are,  respectively,  Poisson’s  ratio  and 
the  modulus  of  rigidity  of  the  conduit.  Shear  and 
pressure  wave  speeds  in  cool  Mount  Hood  andesite 
are  given  by  Murase  and  McBirney  (1973)  as  2700 
and  5500  m/s,  respectively,  and  density  is  given  as 
2500  kg/m3.  The  modulus  of  rigidity  and  Poisson’s 
ratio  of  the  conduit  material  are  then  estimated  as 
1.77  x  10’0  Pa  and  0.34  respectively.  Other  factors, 
not  considered  in  eq  20,  may  influence  the  fluid  wave 
speed.  Jenkner  (1971)  has  shown  that  the  wave  speed 
of  a  fluid  in  a  rectangular  conduit  is  reduced  signifi¬ 
cantly  relative  to  that  in  a  circular  conduit.  Also, 
the  effect  of  trapped  gas  volumes  or  entrained  gas 
bubbles  upon  wave  speed  may  be  important  if  de¬ 
gassing  of  magma  is  suspected  in  low  pressure  regions 
of  the  fluid. 

Initially  the  system  is  pressurized,  but  the  fluid 
is  at  rest.  As  a  result  of  processes  occurring  at  depth, 
additional  fluid  gradually  moves  into  the  system,  and 
the  static  pressure  increases.  Finally  the  increasing 
pressure  causes  stresses  to  develop  in  the  rock  wall 
of  the  conduit  which  exceed  the  strength  of  the  con¬ 
duit,  resulting  in  a  fracture  near  the  closed  end  and 
an  abrupt  outflow  of  fluid.  The  response  of  the 
fluid  in  the  conduit  due  to  this  outflow  is  of  interest 

Neglecting  friction,  the  characteristic  impedance 
of  the  assumed  segment  of  the  Mount  St  Helens 


fluid  system  is  0.0183  using  eq  14  and  20.  The  varia¬ 
tion  of  fluid  system  damping  owith  orifice  impedance 
is  plotted  in  Figure  5.  By  using  Figure  3  to  translate 
system  damping  to  decay  time  of  a  disturbance,  it 
is  clear  that  even  in  the  absence  of  friction,  dissipation 
of  oscillatory  energy  via  transmission  out  of  the  system 
can  be  large. 

A  numerical  solution  of  the  free  vibrational  eigen¬ 
value  problem  was  obtained  to  determine  the  complex 
frequencies  of  the  fluid  system  subject  to  the  frictional 
effects  of  the  high  viscosity  melt.  Circular  frequencies 
were  unchanged  from  the  values  given  in  eq  18  and  19. 
The  lowest  five  harmonics  have  periods  ranging  from 
0.6  to  2.8  s,  consistent  with  our  observations  of  low- 
frequency  oscillations  at  Mount  St  Helens. 

The  magnitude  of  system  damping  o,  including 
the  effects  of  friction,  was  the  same  for  the  first,  third 
and  fifth  harmonics  for  the  single  conduit  fluid  system. 
The  same  was  true  for  the  even  harmonics  of  the  sim¬ 
ple  fluid  system.  The  variation  of  damping  with  ori¬ 
fice  impedance  is  again  given  in  Figure  5.  As  the  ori¬ 
fice  impedance  approaches  the  characteristic  impe¬ 
dance,  the  high  viscosity  and  frictionless  plots  coin¬ 
cide.  This  condition  is  termed  a  matched  line.  In  this 
case,  the  majority  of  the  energy  in  the  oscillation  is 
not  being  reflected  at  the  orifice,  and  the  damping  of 
the  oscillation  due  to  the  transmission  of  energy  from 
the  system  dominates  the  effects  of  friction.  At  the 
highest  and  lowest  flow  extremes,  where  the  ratio  of 
the  impedances  is  much  different  from  one,  most  of 
the  energy  in  the  oscillation  is  being  reflected  by  the 
conduit  termination.  The  effect  of  friction  upon 
system  damping  becomes  apparent,  and  a  lower  limit 
on  damping  is  indicated.  For  a  chosen  value  of  orifice 
impedance,  comparison  of  the  damping  values  given 
by  the  two  curves  indicates  the  importance  of  viscous 
damping  upon  the  overall  damping  of  an  oscillation. 

Figure  5  is  symmetric  about  the  matched  line  con¬ 
dition,  Zc  =  Zor.  Fluid  transients  provide  a  common 
source  mechanism  for  both  pre-eruptive  tremor  and 
tremor  occurring  during  eruptions.  Odd  harmonics  of 
the  conduit,  corresponding  to  a  closed  pipe  condition, 
occur  when  orifice  impedance  is  larger  than  the  char¬ 
acteristic  impedance  of  the  system.  During  pre-eruptive 
activity,  the  conduit  behaves  as  a  closed  pipe.  The 
fluid  outflow  rate  is  relatively  low  and  orifice  impe¬ 
dance  is  large.  Tremors  of  various  durations,  repre¬ 
sented  by  the  right  half  of  Figure  5,  may  accompany 
fluid  movement  The  left  half  of  Figure  S,  where 
orifice  impedance  is  low,  corresponds  to  a  relatively 
large  volumetric  outflow  through  an  open-ended  con¬ 
duit  Even  harmonics  of  the  conduit  occur  during 
these  open-ended,  low  orifice  impedance  conditions. 
Again,  tremors  of  various  durations  can  be  expected. 
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Figure  5.  Variation  of  system  damping  for  the  single  conduit  fluid  system  and  assumed  parameter  values  for  Mount 
St  Helens,  as  a  function  of  orifice  Impedance.  Values  of  o  are  plotted  both  neglecting  and  including  fluid  frictional 
effects.  Small  orifice  Impedance  values  correspond  to  relatively  large  outflow  rates,  while  large  values  reflect  rela¬ 
tively  small  flow  rates.  For  a  matched  line  condition,  Zc=ZOP  damping  Is  very  large. 

The  current  period  of  seismic  activity  at  Mount  durations  as  a  result  of  an  excitation  of  the  near- 

St  Helens  began  in  March  1980.  Since  that  time,  surface  portion  of  the  fluid  system.  The  fluid  dynamic 

episodes  of  volcanic  tremor  with  durations  ranging  source  mechanism  could  be  responsible  for  low-  fre- 

from  several  seconds  to  more  than  an  hour  have  been  quency  vibration  in  any  segment  of  the  volcanic  fluid 

recorded.  As  depicted  in  Figure  5,  the  unsteady  system.  Recent  dome  building  events  and  eruptions, 

flow  source  mechanism  reveals  the  potential  for  a  however,  are  signaled  only  by  near-surface  fluid  tran- 

continuous  variation  of  system  damping  as  a  function  sients  centered  near  the  surface  outlet 

of  the  hydraulic  characteristics  of  the  orifice  or  frac-  The  amplitude  of  volcanic  tremor  at  Mount  St 

ture.  As  orifice  conditions  change,  low-frequency  Helens  has  varied  over  several  orders  of  magnitude, 

earthquake  durations  resulting  from  a  single  distur-  In  general,  tremors  of  smaller  magnitude  have  been 

bance  can  potentially  vary  between  the  shortest  dur-  recorded  since  the  catastrophic  eruption  of  1 8  May 

ation  type  B  event  and  the  lengthy  episodes  of  har-  1980.  The  amplitude  of  a  tremor  is  proportional  to 

monic  tremor  recorded.  Large  differences  in  the  the  amplitude  of  the  pressure  oscillation  in  the  con- 

durations  of  these  signals  mask  the  existence  of  a  duit  Amplitudes  that  vary  over  a  wide  range  are 

common  source  mechanism.  If  outflow  conditions  possible  both  during  and  prior  to  an  eruption, 

are  not  changing  significantly  and  short  duration  Prior  to  an  eruption,  a  limiting  condition  on 

type  B  events  have  recently  been  recorded,  subse-  tremor  amplitude  is  provided  by  the  overburden 

quent  tremor  originating  from  the  same  location  pressure  at  the  point  of  fluid  outflow.  Larger  tremor 

would  likely  be  due  to  either  a  number  of  excitations  amplitudes  are  possible  before  catastropic  failure  if 

of  the  fluid  system  that  are  closely  spaced  in  time  a  large  overburden  is  present  Pressures  greatly  exceed 

or  the  existence  of  some  mechanism  that  is  forcing  static  pressure  during  unsteady  flow  events.  Rock 

the  oscillation.  that  fractures  to  initiate  the  unsteady  flow  is  likely 

The  many  physical  changes  that  have  occurred  to  fracture  further  during  the  transient 

in  the  vicinity  of  the  present  lava  dome  at  Mount  During  the  period  of  23  April-1 7  May  1 980,  the 

St  Helens  potentially  allow  a  wide  variation  of  tremor  north  flank  bulge  at  Mount  St.  Helens  expanded  at  a 


9 


rate  of  approximately  1.5  m  per  day.  As  the  area 
of  the  bulge  was  about  4  km2,  this  expansion  rate 
corresponds  to  an  average  continuous  outflow  of 
approximately  70  m3/s  from  the  conduit  Depending 
upon  orifice  impedance,  oscillatory  flows  of  less  than 
±  100  m3/s  could  produce  initial  pressure  oscillations 
oft  2.4  x  10®  Pa.  Much  larger  pressure  oscillations 
prior  to  1 8  May  are  not  probable  as  the  order  of 
magnitude  of  the  overburden  pressure  of  the  moun¬ 
tain  would  have  been  approached.  Relieved  of  the 
massive  overburden,  static  pressure  and  tremor  ampli¬ 
tudes  are  limited.  Magmatic  eruptions  should  con¬ 
tinue  to  yield  relatively  smaller  volumes  of  material 
until  a  large  overburden  pressure  is  again  developed. 

A  characteristic  of  Mount  St.  Helens  similar  to 
those  of  Hawaiian  volcanoes  is  the  presence  of 
tremors  during  major  volcanic  activity  (Kinoshita 
et  al.  1969).  Tremor  amplitudes  in  both  cases  typ¬ 
ically  subside  as  the  intensity  of  the  activity  dimin¬ 
ishes.  As  an  eruption  wanes  and  mean  flow  dimin¬ 
ishes,  the  magnitude  of  flow  oscillations  can  also 
be  expected  to  diminish.  For  an  orifice  that  is  not 
radically  changing  in  size,  orifice  impedance  is  con¬ 
stant.  As  pressure  and  flow  oscillations  at  the  ori¬ 
fice  are  related  by  the  orifice  impedance  (eq  16),  the 
magnitude  of  pressure  oscillations  at  the  orifice  di¬ 
minishes  with  the  flow.  If  we  assume  that  the  sec¬ 
ond  harmonic  of  the  conduit  is  excited,  the  maximum 
pressure  oscillation  occurs  at  the  mid-length  of  the 
conduit  Its  magnitude  is  related  directly  to  that 
at  the  orifice  by  the  transfer  functions.  It  is  gener¬ 
ally  believed  that  more  than  1  km3  of  material  was 
ejected  from  Mount  St.  Helens  on  18  May  1980. 

The  greatest  portion  of  the  flow  occurred  during  the 
first  9  hours  of  the  eruption.  The  average  flow  rate 
during  that  period  can  then  be  estimated  at  30,000 
m3/s.  Small  pressure  oscillations  at  the  orifice,  cor¬ 
responding  to  a  flow  oscillation  of  ±5,000  m3/s, 
yield  a  pressure  oscillation  of  ±2.4  x  106  Pa  at  the 
pipe  mid-length.  As  the  outflow  oscillations  diminish 
by  an  order  of  magnitude,  the  maximum  pressure 
oscillations  diminish  correspondingly. 

For  an  elastic  pipe,  radial  wall  deformation  in 
response  to  the  oscillatory  pressure  can  be  estimated. 
These  calculations  suggest  that  sufficient  energy  is 
available  via  the  fluid  dynamic  source  mechanism 
to  generate  the  range  of  observed  tremor  magnitudes. 

Preliminary  examination  of  our  seismograph 
records  indicates  a  variation  in  the  spectra  of  both 
type  B  events  and  harmonic  tremor  episodes.  One 
explanation  is  that  different  harmonics  are  excited 
at  different  times.  Other  possibilities  are  that  the 
outflow  location  has  changed,  modifying  the  effective 
conduit  length,  or  that  the  conduit  length  has  in¬ 
creased  due  to  a  hydraulic  fracturing.  Another  pos¬ 


sibility  is  that  different  components  of  the  fluid 
system  are  experiencing  transient  conditions  at  dif¬ 
ferent  times. 

Many  low-frequency  seismic  signals  generated  at 
volcanoes  exhibit  an  emergent  onset  Over  the  first 
few  cycles,  the  magnitude  of  the  oscillation  increases. 
If  this  characteristic  is  source-related,  it  requires  a 
modification  of  the  static  reservoir-conduit-orifice 
fluid  system  component  In  general,  fluid  transients 
generate  the  largest  pressure  oscillations  just  after  the 
disturbance  occurs.  Without  continued  forcing,  these 
oscillations  decay  in  time.  A  diaphragm  inserted  up¬ 
stream  of  the  orifice,  to  account  for  conduit  volume 
change  in  response  to  fluid  pressure  change,  will  pro¬ 
duce  an  emergent  fluid  pressure  oscillation.  An  ori¬ 
fice  thn  is  changing  in  size  with  time  could  also  pro¬ 
duce  an  emergent  pressure  oscillation. 

Beating  is  another  characteristic  observed  in  the 
low-frequency  seismic  output  of  volcanoes  (Fig.  1). 
The  combined  effects  of  different  but  concurrent 
outflows  from  a  segment  of  the  fluid  system  will 
produce  beating  oscillatory  pressures  in  the  fluid. 
Alternatively,  different  segments  of  the  fluid  system 
could  be  experiencing  transient  flow  conditions  con¬ 
currently.  This  concurrent  transient  flow  condition 
occurs,  for  example,  when  two  conduits  of  unequal 
length  are  connected  by  a  crack  which  opens  to  allow 
fluid  passage  when  a  sufficient  pressure  difference 
across  it  develops. 

CONCLUSIONS 

We  have  investigated  the  concept  of  a  fluid  dynamic 
source  mechanism  of  low-frequency  volcanic  earth¬ 
quakes.  The  concept  is  suggested  by  analogous  phys¬ 
ical  processes  occurring  in  three  other  settings:  1)  hy¬ 
draulic  fracturing  of  rock  formations,  2)  low  frequency 
icequakes  in  glaciers  and  3)  auto-resonance  in  hydro- 
power  stations. 

The  analysis  has  linked  long  and  short  duration 
tremors  originating  at  arbitrary  depth,  occurring  both 
prior  to  and  during  volcanic  eruptions,  to  a  common 
source  mechanism.  The  existence  of  tremors  of 
fluid  dynamic  origin  is  independent  of  fluid  proper¬ 
ties  and  conduit  geometry,  although  the  nature  of 
the  tremor  is  affected  by  these  parameters. 

The  fluid  dynamic  source  mechanism  exhibits  all 
of  the  features  generally  attributed  to  the  source 
mechanism  of  volcanic  tremor  including  1)  a  station¬ 
ary  seismic  source,  2)  a  source  location  that  can  vary 
between  tremor  episodes,  3)  a  peaked  frequency 
spectrum  and  4)  a  dominant  frequency  which  can 
change  during  the  course  of  eruptive  activity. 
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The  fluid  dynamic  theory  can  be  applied  in  as 
much  generality  as  is  warranted  by  data  availability. 
We  have  used  features  of  the  recent  eruptive  sequence 
at  Mount  St  Helens  to  demonstrate  the  application 
of  the  linear  theory  of  unsteady  fluid  flow.  With 
more  precise  information  regarding  the  fluid  system 
configuration,  tremor  magnitudes  and  source  loca¬ 
tions,  greater  quantitative  insight  could  be  obtained. 
Even  without  this  information,  the  analysis  has  pro¬ 
vided  several  insights.  Variable  system  damping  can 
result  in  a  broad  range  of  low-frequency  seismic  sig¬ 
nal  durations.  For  a  simple  system,  damping  is 
affected  by  characteristics  of  the  fluid,  conduit  and 
outflow.  It  has  been  demonstrated  that,  for  a  given 
system,  a  change  in  outflow  conditions  is  all  that  is 
required  for  large  changes  in  system  damping  char¬ 
acteristics. 

The  magnitudes  of  tremors  are  related  directly 
to  the  magnitude  of  oscillatory  fluid  pressures.  For 
the  case  of  pre-eruption  tremor,  these  pressures  are 
an  indication  of  the  magnitude  of  the  potential 
eruption.  Tremor  amplitudes  during  an  eruption 
generally  monitor  its  progress,  decreasing  as  the 
eruption  wanes.  Total  fluid  movement  cannot  be 
quantified  from  volcanic  tremor  data  alone  since 
steady-state  or  slowly  changing  fluid  flows  are  aseis- 
mic. 
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